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We investigate the BCS treatment of neutron-proton pairing involving time- 
"gj \ reversed orbits. We conclude that an isospin-symmetric hamiltonian, treated 

pLn \ with the help of the generalized Bogolyubov transformation, fails to describe 

the ground state pairing properties correctly. In order for the np isovector 
pairs to coexist with the like-particle pairs, one has to break the isospin sym- 
metry of the hamiltonian by artificially increasing the strength of np pairing 
^ , interaction above its isospin symmetric value. We conjecture that the np 

T-j- \ isovector pairing represents part (or most) of the congruence energy (Wigner 

^D ' term) in nuclear masses. 
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I. INTRODUCTION 

Pairing correlations are an essential feature of nuclear structure |]I|] . In proton- rich nuclei 
with N K, Z the neutron and proton fermi levels are close to each other and therefore 
the neutron-proton {np) pairing correlations can be expected to play a significant role in 
their structure and decay (for a review of the early work on np pairing theory see Ref. 
0). In contrast, in the heavier nuclei with large neutron excess the neutron-proton pairing 
correlations can be usually neglected. 

There has been a recent revival of interest in the theoretical description of pairing involv- 
ing both neutrons and protons ||^-§[ • This renaissance stems from the advent of experiments 
with radioactive beams, as well as from the application of neutron-proton pairing concepts 
in the description of alpha decay |^ and double beta decay |§,^. However, the theoretical 
treatment is not without a controversy. While intuition and arguments of isospin symme- 
try suggest that the neutron-proton pairing correlations should be as important as the like 
particle pairing correlations in the A^ ~ Z nuclei, the balance between these pairing modes 
is delicate and the standard approximations often fail. 

In order to elucidate what is going on we examine the treatment of neutron-proton 
pairing in the generalized Bogolyubov transformation approach, in particular the role of 
isospin symmetry. The problem at hand is the determination of the ground state of an 
even-even system with the hamiltonian 
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where {jmt) represents the angular momentum, its projection and the isospin projection 



of the single-particle (s.p.) state created (annihilated) by the operator a] 
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The three coupling constants Gu' (we assume that Gu' 



Gt't) characterize the monopole pairing interaction. The interaction couples only states in 
time reversed orbits, but allows an arbitrary combination of the isospin projection indeces. 
Obviously, when isospin symmetry is imposed, the s.p. energies become independent of the 
isospin label t, and Gnn = Gpp = Gpn = Gnp = G. The interaction then describes the 
isovector T = 1 pairing. However, as will be seen below, it is advantageous to keep the 
general form of the hamiltonian (|I|) . 

One can find the exact ground state of (|I]) in the simple case of one or two level systems. 
However, in the general case of a multilevel system the dimension increases exponentially 
and therefore the standard procedure is to use the generalized Bogolyubov transformation 
approach in the form |jTO|, where the quasiparticle operators are related to the particle 
operators by 
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Here j denotes the full set of quantum numbers of a s.p. orbit, and the indeces "1" and "2" 
are the quasiparticle analogs of p or ra, i.e. of the corresponding isospin projections. The 



transformation amplitudes Uikj and Vikj with i ^ k describe the neutron-proton pairing. 
They are, in general, complex. We refer to [^ and [jlO| for the unitarity conditions which 
Uik^j and Vikj have to obey, as well as for the relation between the amplitudes and the gap 
parameters Ap, A„, and A„p. 

To find the ground state we minimize the quantity Hq, the expectation value of the 
hamiltonian in the quasiparticle vacuum, while simultaneously obeying the unitarity condi- 
tions and the usual conservation (on average) of the number of neutrons and protons. (This 
procedure is equivalent to demanding that the "dangerous graph" term if 20? which creates 
or annihilates a pair of quasiparticles, vanishes.) We use the Newton-Raphson method |TT 



and check, by comparing to the "standard BCS" solution for Gnp = 0, that the ground state 
energy is lower than in the state without the neutron-proton pairing. The procedure allows 
us to find at the same time the gain in the ground state binding energy associated with the 
neutron-proton pairing. 

II. ISOSPIN SYMMETRIC HAMILTONIAN 

The hamiltonian (|1]) with ejp = ej„ and Gnn = Gpp = Gnp describes the isovector 
pairing, in which all three kinds of pairs {nn, pp, and np with T = 1) are treated equally 
on the interaction level. One expects then that in an even-even nucleus with N = Z the 
corresponding gap parameters should be the same for all three possible pairs. 

In fact, in the exactly solvable manifestation of this hamiltonian, in which there is only 
one s.p. state of degeneracy 2Q, this is indeed the case. Defining the pair creation operator 
as 

^t^ = z^ y^jm^j^t^ ? (3) 

where tz is the corresponding isospin projection, the quantity related to the pairing gap Aj^ 
is the ground state expectation value Mt^ = (Sl^St^). (We calculate the "gap" Aj^ from the 
expression A/t^ = A^^/G^^ valid up to the terms 1/Q. This relation, however, fails for full 
shells.) As shown in 0, based on the earlier work on this 5*0(5) model, one can obtain 
analytic expressions for A/j^. Indeed, when N = Z and both are even, all three values of 
A/t^ are equal, and when N — Z increases the A/q, and therefore also A„p, sharply decreases, 
while the other two A/t^=±i remain the same or increase with N — Z. We expect that this 
behavior is "generic", i.e. survives even in the case of more than one single particle level. 

The conjecture above is supported by the analysis of the two level model |T2|. If we 
restrict ourselves only to the states with seniority zero it is easy to construct the corre- 
sponding hamiltonian matrix which has then very manageable dimensions even for large 
Q. For completeness we give the expressions for the corresponding matrix elements in the 
Appendix. 

Unlike the exact solutions described above, the generalized Bogolyubov transformation 
approach gives a very different results in the isospin symmetric case. It has been known for 



some time |13] that in the one- level case the approach leads to no np pairing when N > Z. 
Our calculation shows that this is true generally. There is, in fact, a no mixing theorem: The 
ground states of even-even nuclei with N > Z have vanishing A„p when isospin symmetric 
hamiltonian is used. For N = Z nuclei there is still no mixing. But in that case there are 



two degenerate minima of the energy; one with nonvanishing A„ = Ap and A„p = 0, and 
the other one with A„ = Ap = and A„p 7^ 0. (This conclusion was also reached in 0] for 
the one-level model, and in in the more general case.) 

We see, therefore, that the generalized Bogolyubov transformation fails to describe cor- 
rectly the treatment of isovector pairing correlations. However, since we expect, as stated 
above, that the effect of np pairing decreases fast with increasing N — Z, the standard BCS 
theory is still applicable for most heavier nuclei where iV — Z is relatively large. 

III. BREAKING THE ISOSPIN SYMMETRY 

Let us consider now what happens when the requirement of isospin symmetry is relaxed, 
i.e., in the hamiltonian (|I]) one allows different coupling constants Gnn 7^ Gpp 7^ Gnp, and 
possibly different single particle energies for neutrons and protons. It was shown already 
thirty years ago |]TD| that such a hamiltonian, treated using the generalized Bogolyubov 
transformation (j^), results in nonvanishing Anp in nuclei with N c::^ Z. 

Hamiltonian which breaks isospin symmetry leads, naturally, to eigenstates that do not 
have a definite value of isospin. Since the quasiparticle vacuum mixes states with different 
particle number, and therefore also with different isospin, even for an isospin conserving 
hamiltonian, it is perhaps worthwhile to explore effects associated with such a more general 
situation. 

It is straightforward to treat the hamiltonian (|1]) exactly in the one or two level models; 
the corresponding hamiltonian matrix can be calculated using the formuae in |T^ (see also 
Appendix). The corresponding eigenstates are no longer characterized by isospin T. Instead, 
all isospin values between T^ = T^in = [N — Z)/2 and T^ax = {N + Z)/2 contribute to the 
wave function. The ground state energy Egg of a one level system with Gnp 7^ Gpair = Gnn = 
Gpp decreases monotonically with increasing Gnp- However, the binding energy gain between 
a system with no neutron-proton interaction (and therefore A„p = 0), and the system with 
pure isovector interaction (and A„p 7^ 0) is only of the order 1/^2, 

AE = EgsiGnp/Gpair = 1) ~ Egs{Gnp = 0) = —GpairZ/2 , (4) 

compared to the leading term —Gpair^{N + Z)/2. Moreover, the exact wave function of the 
ground state corresponding to the isovector pairing hamiltonian with G^p/ Gpair = 1 can be 
obtained from the ground state of the isospin violating hamiltonian with any Gnp/ Gpair 7^ 1 
by simply projecting onto a state with isospin T = Tmin- This is an exact statement which 
follows from the uniqueness of the zero seniority state with given (A^, Z,T). 

In Fig. 1 we show the exact and BCS gap parameters for the one level system as a 
function of the ratio Gnp/ Gpair- The degeneracy Q and the pairing strength Gpair are chosen 
in such a way that they resemble the situation in finite nuclei discussed later on. One can 
see that the two methods give qualitatively similar results. They agree with each other 
quite well, with the exception of the narrow region near the "critical point" of the BCS 
method (for the plotted case this point is at Gnp/ Gpair = 1-05). As usual, the BCS method 
is characterized by the sharp phase transition while the exact method goes more smoothly 
through the "critical point" , as it must for a finite system. Nevertheless, the basic similarity 
is apparent. 



It is now clear that the failure of the BCS method to describe the neutron-proton pairing 
in the isospin symmetric case is not a fundamental one. It is related to the abrupt phase 
transition inherent in the BCS. The isospin symmetric value Gnp = Gpair is less than the 
critical value needed for the phase transition from the pure like-particle pairing to the sit- 
uation where both like-particle and neutron-proton pairs coexist. Since, as stated earlier, 
the quasiparticle vacuum breaks isospin anyway, it should not matter much that the isospin 
violation is also imposed on the hamiltonian level. We have to choose, however, the proper 
value of the coupling constant Gnp- 

The natural way to fix the ratio Gnp/Gpair is in nuclei with N c::^ Z where one can use 
the arguments of isospin symmetry to estimate the gap A„p. We show in Fig. 2 again 
the comparison between the exact and BCS gaps in the degenerate case, now as a function 
oi N — Z. There are basic similarities between the two situations, but the quantity Anp 
decreases more rapidly with N — Z in the BCS case than in the exact case. We believe that 
this feature is related to the approximation involved in relating the gap Aj^ to the ground 
state expectation value Aft^ in the exact case. What is clearly visible in both cases, and 
intuitively obvious, is the tendency of the A„p to decrease with increasing N — Z. This 
tendency have been noted many times before, see e.g. 0J^. In particular Ref. [Q have 
shown that in the BCS approach for real nuclei, and with the ratio Gnp/Gpair fixed so that 
at A^ ~ Z the gap A„p has reasonable value, the effect of neutron-proton pairing disappears 
at N - Z > 6. 



IV. NEUTRON-PROTON PAIRING AND CONGRUENCE ENERGY 



It has been known for some time, and reiterated recently in Refs. [pTStplBI that nuclear 
masses can be well approximated by a smooth Thomas- Fermi formula, shell correction terms, 
even-odd terms (i.e. standard like-particle pairing), and the congruence energy (or "Wigner 
term") which depends on |iV — Z\/A. The least squares fit to the congruence energy gives 

m 



C{I) = -10e-^-2|^l MeV; I = {N - Z)/A . (5) 

Satula and Wyss suggested that the neutron-proton pairing interaction in the isospin 
T = (and therefore J 7^ 0) channel provides a microscopic origin of the Wigner term. In 
their work [Q there is no neutron-proton T = 1 pairing because they use isospin symmetric 
hamiltonian and thus are missing the neutron-proton T = 1, J = pairing effects as ex- 
plained above. For the microscopic mechanism of the Wigner term in Ref. [^ it is necessary 
that the strength of the T = pairing (called x^^^ there) be above the critical value x^^^ 
= 1. 

Here we consider only pairing interaction in time-reversed orbits and thus the possible 
T = pairing is not present. However, with the treatment of the T = 1 neutron-proton 
pairing suggested here we have an alternative (or additional) microscopic contribution to 
the congruence energy, caused simply by the neutron-proton T = 1 pairing. To see how 
this alternative would work we compare in Fig. 3 the extra binding energies associated with 
the np T = 1 pairing with the congruence energy C{I) above. (Since we are interested 
mainly in the \N — Z\j A dependence we shifted the C(/) curve by 7 MeV.) As in Ref. H 



the extra binding energy associated with pairing decreases with |/| faster than the empirical 
congruence energy C(|/|), at least for the degree of isospin symmetry breaking we have 
chosen. Obviously, the quantity 

5Bnp = EiGnp = l.25Gpa^r) ' E{Gnp = 0) (6) 

shown in Fig. 3 has qualitatively similar dependence on (N—Z) as the congruence energy, eq. 
(|). (The single particle level scheme end energies of Ref. [^ were used in our calculation.) 
It remains to be seen, however, whether this contribution to the congruence energy is the 
only one (or at least the most important one). Also, further investigation of the dependence 
of the gap A„p on the neutron excess is needed. It is possible that the BCS treatment with 
the fixed ratio Gnp/Gpair exaggerates the slope of that dependence. 

V. CONCLUSIONS 

We have shown that treating pairing properties of a system of interacting protons and 
neutrons with the help of the generalized Bogolyubov transformation requires special care. 
The corresponding system of equations allows, in principle, three different solutions. There is 
the trivial "normal" solution with no pairing whatsoever. But there are also two competing 
solutions with pairing correlations present. One, in which there are no neutron-proton 
pairs, corresponds to a product state with the neutron-neutron and proton-proton pairs not 
communicating with each other. The other solution corresponds to a system in which the 
like-particle and neutron-proton isovector pairs coexist. When the generalized Bogolyubov 
transformation is used, there is a sharp phase transition between these two paired regimes, 
with the critical pairing strength Gnp/Gpair somewhat larger than unity. 

Thus, if one wants to describe the neutron-proton pairing using the quasiparticle trans- 
formation method, one has to break isospin symmetry at the hamiltonian level. We propose 
to fix the unknown degree of isospin breaking in such a way that the gap A„p in iV ~ Z 
nuclei is reasonable, i.e. comparable to the gaps A„„ and App. With this assignment all 
traces of the isovector neutron-proton pairing disappear for N ~ Z > 6 in real nuclei. 

We conjecture that the isovector neutron-proton pairing represents one of the main con- 
tributions to the congruence energy (Wigner term) which describes the additional binding 
of nuclei with {N - Z)/A ~ 0. 
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APPENDIX A: 

For the exact diagonalization in the space of pairs with zero seniority we use the basis 
I A/", T, Tj), where A/" is the total number of pairs, T is the isospin and T^ its projection. The 



necessary matrix elements of the pair creation and annihilation operators are 

(Ar+l,T+l,T, + t,|4jAr,T,T,) = 

■[(T + l)(2f] - AT - T)(T + AT + 3)/2]i/^ (Al) 



(TTJt.|T + l,T. + t,)^,^ , .^.on ^^ r^^^rr. , .. , o^ /01I/2 



V2T + 3 



-[T(2f] + l-Ar + T)(Ar-T + 2)/2]i/^ (A2) 



(TTJt^lT- l,Tj + t2)^^^^^ ,1 A/- , T-\/A/- T- , 0\/0ll/2 



V2T-1 

{Af-l,T+l,T,-t,\StJM,T,T,) = 

i^^^i^^^^i|^^^^[(T + l)(2fi + 3-Ar + T)(Ar-T)/2]V^ (A3) 



■[T{2Q + 2-Af-T){Af + T+l)/2Y^^. (A4) 



(TT^lt^lT 1,T2 tz) frr^cn, ,0 A/- T^\/A/- , T^ , 1\/Oll/2 



V2T- 1 

For the case of two levels the basis is |A/i, A/2, Ti, T2, T12, T22) with the obvious constraint 
A/i + A/2 = A/" (total number of pairs) and a similar one for T^. The hamiltonian matrix 
is easily constructed from the expressions above, with terms diagonal in A/i,A/2 and terms 
where A/i — > A/i ± 1, and A/2 — > A/2 =F 1. In addition there is a diagonal shift of 2eA/2, where 
e is the splitting of the single particle states. 
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FIGURES 

FIG. 1. The pairing gaps for the one-level case with 0, = 11, A^ = 6 and Z = 4. Gpair = 0.242 
MeV was used and the results are plotted as a function of the ratio Gnp/Gpair- The upper panel 
is for the exact solution with gaps determined as described in the text. The lower panel is for the 
BCS solution. 

FIG. 2. The pairing gaps for the one-level case with fi = 11, as function of the neutron excess 
N — Z. Both panels are for Z = 4 while the neutron number is varied; Gpair = 16/ (A^ + Z + 56) 
MeV. The exact solutions for Gnp/Gpair = 1-1 are in the upper panel. The BCS solutions for 
Gnp/Gpair = 1-25 are in the lower panel. 

FIG. 3. The extra binding energy associated with the neutron-proton pairing, eq.(|6|) is depicted 
by the dashed lines for the series of Ge (upper panel), Se (middle panel) and Kr (lower panel) 
isotopes. The full line, shown for comparison, represents the congruence energy, eq.(^), shifted by 
7 MeV. 



> 



1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 \— 
1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 



1 1 1 1 1 1^' 




/ 
/ 
/ 


"^ \ 


^ 


/\ - 


^N^ N^ 


/^^\ - 


/ 
/ 

y 


— -^^F \ \ \ \ H 


_ 


y 


~ 


1 1 


1 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

np/ pair 



> 




N-Z 



> 






Oh 




-6-4-2 2 4 6 



N-Z 



